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Routers, switches, and repeaters are essential 
components of modern information-processing 
systems. Similar devices will be needed in fu- 
ture superconducting quantum computers. In 
this work we investigate experimentally the time 
evolution of Autler-Townes splitting in a super- 
conducting phase qubit under the application of 
a control tone resonantly coupled to the second 
transition. A three-level model that includes 
independently determined parameters for relax- 
ation and dephasing gives excellent agreement 
with the experiment. The results demonstrate 
that the qubit can be used as a ON/OFF switch 
with 100 ns operating time-scale for the reflec- 
tion/transmission of photons coming from an ap- 
plied probe microwave tone. The ON state is re- 
alized when the control tone is sufficiently strong 
to generate an Autler-Townes doublet, suppress- 
ing the absorption of the probe tone photons and 
resulting in a maximum of transmission. 

In the optical regime, a variety of effects that can be 
employed for switching have been demonstrated, for ex- 
ample photon blockade by single atoms ^, parametric in- 
stabilities in nonlinear optical wave-mixing^, heralded 
single-photon absorption by one trapped ion^, single- 
atom dynamic control of light in microresonator pho- 
ton turnstiles^, and population inversion in single dye 
molecules^. Optical switching using the phenomenon 
of electromagnetically induced transparency has been 
realized in gases of magneto-optically trapped ^^Rb ; 
recently, it has become possible to use for this pur- 
pose only a single ^^Rb atom in a high- finesse cavity . 
With improved dissipation and decoherence, it could be 
possible to use electromagnetic-induced transparency ef- 
fects to build similar devices with quantum supercon- 
ducting circuits'^ Other applications employing multi- 
level systems - for example state preparation, emulation 
of large quantum- number spins , single-shot fast quan- 
tum nondemolition readout techniques'^, implementa- 
tion of control Z-tt gates in certain quantum computing 
architectures , single-qubit microwave amplifiers , etc. 
- can be envisioned. 

In this paper we demonstrate the dynamical operation 
of a phase qubit as a microwave photon-absorbing device 
that can be switched ON or OFF in approximately 100 ns 
by using an external control field. The functioning of our 
device is based on an effect closely related to electromag- 
netically induced transparency: under the application of 
a relatively intense continuous field (the control field) to 



the second transition of the qubit, the spectral line of 
the first transition splits into a doublet. This is called 
Autler-Townes effect'^, and it has been recently observed 
in superconducting quantum systems as well'^ Here 
we operate our device in time-domain, which allows us to 
observe transient effects that occur until a steady state 
is reached. We observe how the doublet forms when the 
control field is suddenly applied, and how the spectral 
line of the first transition appears back from the merg- 
ing of the Autler-Townes peaks when the coupling field 
is switched off. 

Besides the fundamental importance for understanding 
the dynamics of the Autler-Townes effect, our experiment 
has potential applications in the field of microwave pho- 
tonics, for example for fast quantum switches that can 
be integrated in superconducting quantum-processor ar- 
chitectures. This can be achived by realizing that the 
"control" microwave tone applied to the qubit deter- 
mines the absorption rate for a "probe" microwave tone. 
Indeed, when the control field is OFF, the third level 
does not play any role and the device absorbs the inci- 
dent probe field radiation resonant to the lowest transi- 
tion, and its first excited level is populated. If such a 
two- level system is embedded into a transmission lines, 
it has been predicted theoretically in quantum optics'^ 
and demonstrated experimentally that the qubit will 
reflect back the incoming photons by resonant fluores- 
cence, provided that the intensity of the probe field is 
not too high (this saturates the qubit). When the control 
field is ON, the absorption from the qubit is suppressed 
almost completely - the system becomes "transparent" 
to the incident probe radiation and the transmission is 
maximal. We calculate theoretically the expected reflec- 
tion coefficient both in the transient and in the steady- 
state regimes, and we also analyze the performance of 
the switch at different strengths of the control and probe 
tones. Using our theoretical model, we are able to present 
a set of general quidelines for the design of supercon- 
ducting switches based on the Autler-Townes effect - we 
show that high anharmonicity, low probe powers, and low 
dephasing rates lead to an effective suppressions of the 
reflection. 



Results 

Our phase qubit^^ consists of an r.f. SQUID with loop 
inductance L, junction capacitance C, and Josephson en- 
ergy which is read out by another d.c. SQUID fabri- 



2 




FIG. 1: Image of the Josephson phase qubit sample 
used in the experiment. A Al/AlO^/Al junction with an 
area of ~ 14 /xm^ was patterned on a sapphire substrate by 
standard optical lithography. SiN^ was used as an insulator 
between different metallic layers. 
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FIG. 2: Schematic of the pulse sequence for the dy- 
namic Autler-Townes effect. The qubit is probed by a 
continuous microwave field (uOp^Qp), with the frequency Up 
being swept around the first frequency transition cjiq. Each 
measurement sequence starts at t = ns with the system in 
the stationary state as given by the application of the probe 
field. At t = 150 ns we apply the coupling pulse (cjc^c), 
exactly resonant with the second transition, ujc = 0J2i- 



cated on-chip in the proximity of the qubit and coupled 
inductively to it (see Fig. 1). The qubit can be biased by 
an externally applied magnetic flux and can be coupled to 
other circuit elements such as transmission lines and res- 
onators. The device forms a multilevel quantum system 
with the first three levels denoted by |0), |1), and |2). 
These energy levels can be addressed using microwave 
fields closely resonant to the first |0) |1) and the sec- 
ond |1) 1 2) transitions (the |0) |2) transition is 
almost forbidden). 

A detailed description of the functioning of the device 
has been given elsewhere^ ^. The Hamiltonian of the sys- 
tem is 
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where Q is the charge on the capacitor formed by the 
junction, <l> is the flux variable, <l>o = h/2e = 2.067 x 
10~^^ Wb is the magnetic flux quantum, and <l>ext = 
^dc + ^rf (^) is the total magnetic-flux component of the 
externally applied d.c. and r.f. fields. In this system, 
the r.f. tones applied have angular frequencies ujp (probe 
tone) and uJc (control tone) close to the two resonant fre- 



quencies ujio of |0) ^ |1) transition and UJ21 of |1) ^ |2) 
transition, respectively (see Fig. 2). The first two tran- 
sitions were determined in independent single-tone and 
two-tone spectroscopy measurements^^, and they had the 
values cjio = 27r x 8.135 GHz and UJ21 = 27r x 7.975 GHz. 
The probe and control fields are supplied to the device 
through an on-chip coplanar waveguide transmission line. 
The Rabi frequencies of these fields are denoted by ftp 
and respectively, and they are directly proportional 
to the corresponding probe and control field amplitudes. 

Experimental data. Here we demonstrate that this 
three-level superconducting device can be operated dy- 
namically, in time domain. In Fig. 2 we present the pulse 
sequence applied in order to use the system as a switch 
based on the Autler-Townes effect. In Fig. 3 we show 
the results of the experiment. The ON state is charac- 
terized by the suppression of excitations to the state |1) 
normally caused by absorbtion of photons from the probe 
field. It is produced by a large value of the control field 
amplitude (effective Rabi frequency ftc)- The state OFF 
corresponds to a low enough value of the control field so 
that excitations to the state |1) are allowed. Although 
we do not measure the transmission (reflection) of pho- 
tons directly, the data demonstrates unambiguously that 
the device can be used to modulate the transfer of pho- 
tons by controlling the absorbtion rate of the qubit. The 
measured state occupation numbers show clearly the dy- 
namical process of formation of the Autler-Townes dou- 
blet, how it reaches the stationary state, and how the two 
peaks finally collapse into a single spectroscopic signal at 
the first transition frequency when the coupling field is 
switched off. In both the on-set and the switch-off stages 
of the effect, the timescale for the system to reach the 
steady state is of the order of 100 ns, in agreement with 
the decoherence times in our system (see the theoretical 
model below). 

Numerical simulation. In order to simulate the 
experimental results including the effects of decoher- 
ence, we use the standard Liouville superoperators in 
the Markov approximation. For our system, described 
by the Hamiltonian equation (1) and placed in a reser- 
voir of electromagnetic modes at zero temperature (see^^ 
for a detailed derivation) , the Markovian master equation 
for the density matrix p in the rotating frame reads 



(2) 



The effective Hamiltonian H is obtained by truncating 
equation (1) to the lowest three levels and performing a 
rotating wave approximation. It can be put in the form 
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^{n, + lA5npe'^'^'^)ai2] + h.c. (3) 

Here we have taken an external r.f. magnetic field 
of the form ^rf(^) = ^pCOs{ujpt + ^) + ^cCOs(a;c^). 




8.05 



FIG. 3: On-set and extinction of the Autler-Townes 
effect. The figure shows the experiment ahy determined sum 
of excited-states level-occupancy probabilities, P|i) +P\2)^ ob- 
tained by applying the pulse sequence from Fig. 2. The arrow 
indicates the frequency of the first transition, ujiq/2ti = 8.135 
GHz. 

This produces an effective coupling to the qubit 
^(^) = gpCos{ujpt + (/))+ gcCos{uJct)^ where 
9c = -{^c/L)y^l/ChuJo and gp = -{%/ L)y^l/ChuJo 
(where ujq is the Josephson plasma frequency) resulting 
in Rabi frequencies ^}p = 0.69gp and = 9c (see 
Ref. 21 for a detailed derivation). Here we have used 

the operators aij ^= (which generalize the Pauli 

matrices to a many- level system), and we have defined 
the detunings Ap = ujiq — uOp^ Ac = 0021 — cjc, and 
5 = ujp — uOc The origin of the qubit-field coupling term 
(the second line in the Hamiltonian equation (3)) can 
be understood easier if one treats the device as a 
harmonic oscillator interacting with fields; in this 
case, the corresponding terms in the Hamiltonian read 
{h/2) [{l/V2){gpe'^ ^gce-'^')cjQi + (^c + ^pe^^'+^^)ai2] + 
h.c. The result of including anharmonic corrections is 
that from the previous expression is replaced by 1.45 
and 1/V^by 0.69. 

The relaxation part is given by^^ 

-F^^ (2cri2pcr2i - CF22P - pcr22) 

-^n (e"^^Voipcr2i + e'^Vi2pcrio) , (4) 
and the dephasing is 

>Cdep[p] = (2criipcrii - (Jup - pan) 

+ ^ (2cr22pcr22 - Cr22P - ^^22) • (5) 

Here p is the density matrix in a doubly-rotating frame"^^, 
the interlevel relaxation rates between |1) |0) and 
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Occupation probabilities P,.. + P|2.(simulation with phase-averaging) 




FIG. 4: Simulated occupation probabilities for the sum 
of the occupation probabilities for the states |1) and 
1 2) with averaging over the phases 0. This quantity is 
measured directly in the experiment (see Fig. 3). 



1 2) |1) are denoted as Tio and r2i, respectively, 7^0, 
72^0 are intrinsic (pure) dephasing rates of the states |1) 
and 1 2), and n = vTToI^- The dissipation parame- 
ters were either measured directly in independent exper- 
iments {e.g. by exciting the qubit and measuring the de- 
cay time) or extracted from spectroscopy data^^'^^, and 
for this sample they were r2i = 27rx 11 MHz, Tio =27rx7 
MHz, 7fo = 27r X 7 MHz, 71^0 = 27r x 16 MHz. 

In Fig. 4 we present the result of the numerical simu- 
lation using equation (2). For this simulation, the ampli- 
tude of the probe field corresponds to Rabi frequencies 
Q.p = 271 X 3.45 MHz; for the control field, the corre- 
sponding values were = 27r x 60 MHz in the ON state, 
and a much smaller value = 27r x 5 MHz in the OFF 
state. The calibration of the probe and control field am- 
plitudes was done in independent measurements of the 
Rabi frequency between the levels |0) and |1), and be- 
tween the levels |1) and |2) respectively. In the absence 
of the control field, the level |2) is empty while the level 
|1) is populated when the probe field is resonant with 
the |0) |1) transition. After the control field has been 
switched on, the original high occupation peak centered 
at C(Jio splits into the Autler-Townes doublet. These ob- 
servations are in good agreement with the experiment 
(Fig. 3). The simulation also captures two additional 
features of the experimental data: the asymmetry in the 
Autler-Townes peaks and the appearance of a plateau 
with nonzero occupation probabilities even when ujp is 
far off-resonant from ujiq in the ON state. Both of these 
effects result from the existence of cross-coupling terms in 
the Hamiltonian equation (3), see the Discussion below. 

In the experiment, the phase difference (j) between the 
probe and the control field cannot be maintained at the 
same value for each realization of the sequence presented 
in Fig. 2. To model this, we calculate the probabilities 
P\i) ^ P\2) foi" several values of (j) randomly distributed 
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in the interval [0, 27r] and we average over the results. 
More considerations about the effect of the cross-coupling 
terms are delegated to the Discussion. 

We note that the power levels used in this experiment 
correspond to operation close to the level of single mi- 
crowave photons. To estimate the average number of 
photons involved, consider a generic process (without 
dissipation) of resonant Rabi oscillations with Rabi fre- 
quency r^R. The probability for a two- level system to 
be excited in a time r is sin^(l]Rr/2). Therefore, after 
a time tr = tt/I^r the system is with certainty in the 
first excited state, meaning that it has absorbed exactly 
one photon. The corresponding power is hw/rp. In our 
case, the power levels correspond to a Rabi frequency 
Qp = 27r X 3.5 MHz, resulting in Tp = 143 ns (and a 
power of hwio/Tp = 38 aW = -134 dBm), and for the 
control field in the ON state we have = 27r x 60 MHz, 
yielding Tc = 8.3 ns (and a power of huj2i/rc = 0.6 fW 
= -122 dBm). Thus, during the time between t = 150 
ns and t = 300 ns there were on average about 1 pho- 
ton present in the probe beam and about 18 photons in 
the control beam. Moreover, as shown in the next sec- 
tion, increasing the input power level saturates the qubit 
and switching does not occur anymore. These low power 
levels are typical for measurements of superconducting 
qubit s, therefore the devices based on the mechanism de- 
scribed in this paper can be seamlessly integrated in any 
circuit QED architecture. 



Discussion 

The effect of the cross couphng terms. In the 

effective Hamiltonian equation (3), there are oscillating 
terms with factors ex.p{±iSt)^ which result from cross 
couplings of the coupling field into the |0) |1) tran- 
sition and of the probe field into the |1) |2) transition 
(see also the detailed derivation in Sec. H and Sec. HI of 
Ref. 21). In the real experiment, it was not possible to 
ensure that the phase between the coupling and the probe 
fields remained the same for every run of the measure- 
ment sequence. To obtain Fig. 4, we have averaged over 
0, which is what occurred effectively in the experiment as 
well. It is however instructive to see what would happen 
for a constant phase (or a single-shot measurement), and 
for the case in which there is no cross-coupling at all. 

In Fig. 5(a) we show the result of simulations with 
the same master equation as equation (2) but with the 
relative phase between the probe and the coupling fields 
fixed at = 0. One can see the appearance of fringe-like 
structures in the spectrum. This figure can be regarded 
as a snapshot corresponding to a single run of the pulse 
sequence. The plateau seen in the previous spectra in the 
region where the coupling field is on and the probe is far 
off-resonant still exists; also the asymmetry between the 
Autler-Townes peaks is preserved. At different (^'s, these 
fringes appear in different places; thus when averaging 
over (j) they will be washed away. 



To get even more understanding of the effects of the 
coupling terms, let us analyze what happens if there is 
no cross-coupling at all. This would be the case for a 
strongly anharmonic system. In this case we can approx- 
imate the effective Hamiltonian as 

= ^ [Apau + (A^ + Ac)a22] + ^ [%ctoi + ^^c^i2]+h.c. 

(6) 

This Hamiltonian is widely used in quantum optics for 
studying electromagnetically induced transparency and 
coherent population trapping . 

Fig. 5(b) shows a time-domain simulation similar to 
that in Fig. 5(a) but with the Hamiltonian H' . The 
time-oscillating terms in £rei[p] are similarly neglected. 
Several features can be noticed immediately. The Autler- 
Townes peaks are this time symmetric, the occupation 
probabilities at ujp far away from Autler-Townes split- 
tings are smaller, and we don't have a plateau anymore 
in the ON state (from 150 ns to 300 ns). The latter two 
features can be explained intuitively in a simple manner: 
the coupling field tends to introduce an additional pop- 
ulation of level |1) when it couples into the |0) |1) 
transition. This effect does not depend on the frequency 
of the probe field, which is why it is appears as a plateau. 

Operation as a quantum switch. 

The capability of distributing entangled microwave 
photons in superconducting quantum networks is 
an essential ingredient for building future quantum 
processors'^. An immediate application of the effect pre- 
sented above is to realize a quantum switch for controlled 
routing of microwave photons. These devices can be re- 
alized by embedding the qubit into an open transmission 
line and using a vector network analyzer to measure the 
scattering coefficients of the probe field, yielding the cor- 
responding reflection and transmission coefficients. In 
the stationary regime, these measurements have been 
done for flux qubits and transmons''^. Here we discuss 
the general dynamic characteristics of these devices, with 
exemplification by using the parameters corresponding to 
our qubit. 

For qubits embedded into one-dimensional transmis- 
sion lines, the reflection coefficient r of the probe tone is 
proportional to the off-diagonal element pio = Tr[pcrio]^: 

^rio 

r = ro—pio, (7) 

where tq is the maximum reflection amplitude^^. Note 
that two conventions are used in the literature, one that 
considers the reflection as positive and the other in which 
it is taken negative. Our equations do not depend on the 
convention used, but we assume that the same convention 
is applied both to r and to tq. In the case of our exper- 
iment, the switching occurs between two steady states, 
and Popp, which are reached after ^100 ns of tran- 
sients. As a figure of merit for switches, the OFF/ON 
power ratio between the reflected powers in the OFF 
state and that in the ON state Rqff/Ron must be as 
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FIG. 5: Simulated P|i) + P\2) with = 0, and (a) with cross-coupling, (b) without cross-coupling. 



large as possible. In the simulations below we will use pa- 
rameters corresponding to our experiment; however, we 
stress that the theoretical results are generic and can be 
applied to any type of superconducting quantum multi- 
level system. 

Cross- coupling terms and anharmonicity We first cal- 
culate numerically the off-diagonal density matrix ele- 
ment p[^Q^ and the reflection coefficient obtained from 
equation 7 in the steady state. For these calculations 
we take both the probe and control fields resonant with 
the respective transitions {Ap = and Ac = 0). This 
corresponds to the experimental data along the arrow 
from Fig. 3. In Fig. 6(a) and (b) we present the normal- 
ized reflected power R = |r/rop simulated with cross- 
coupling (using the Hamiltonian iJ), and respectively 
without cross-coupling (using the Hamiltonian H'). The 
effect of the cross-coupling terms become apparent es- 
pecially at high values of ftc] as this field couples into 
the |0) — |1) transition, the reflected power is not turned 
down completely. We have checked numerically that this 
is due to the real part of p^Q^^ acquiring a non-zero value, 

(st) 

while the imaginary parts of pi^ ^ with and without cross- 
coupling do not differ much. For example, without cross- where 710 = F 



Eq. (3), if S is large, the cross-coupling terms will oscil- 
late fast and average to zero). Some types of qubits, for 
example the flux qubit, do have a high degree of anhar- 
monicity by design; for phase qubits there exist as well 
proposals for increasing the anharmonicity . 

Saturation If the probe power is too large, from Fig. 6 
we can see that, in the OFF state, the power reflection R 
is significantly reduced. This is because in the OFF state 
the system is essentially a two-level system exposed only 
to the probe field radiation. In this case, if Qp is much 
larger than 710, the stationary state of the system in the 
{|0), |1) subspace is maximally mixed, and the reflection 
coefficient drops to zero. This shows that increasing the 
power of the probe field would not improve the charac- 
teristics of a switch. 

Dephasing In the absence of cross-coupling^s and for 
small enough values of the probe field, (1^^/710) <C 1 and 
(Qp/Qc)'^ <C 1, analytical results for for the off-diagonal 
element p^Q*^ in the steady-state are available^^. 
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state, and Re[p^Q^^] = 0, lm[p[^Q-'^ 



with cross-coupling, Re[p^Q^^] r 
in the OFF state, and Re[p['^^^ 



-0.01, lm[p[f] 
-0.1,Im[p^'^^i 



-0.19 in the OFF 
-0.02 in the ON 

state, which give a power ratio Rqff/Ron of 90; whereas 

-0.19 

^10 J -L^^HA^lO J ^ -0.04 

in the ON state, which give a ratio Rqff/Ron of only 3. 

This means that, in the state ON, the normalized 
power reflection R = |r/rop with cross-coupling is al- 
ways larger than that without cross-coupling, and the 
additional power reflection due to the cross-coupling re- 
duces the on/off ratio, Rqff/Ron^ of the switch. We 
conclude that, in order to have an efficient switch, it is 
advantageous to increase the anharmonicity of the qubit, 
which would reduce the cross-couplings (for example in 



27r X 14 MHz and 720 = 



720 = 27r X 27 MHz. In Fig. 7, the real and imaginary 
parts of the reflection coefficient obtained with pi^^ of 
equation (8) are plotted for different values of Qc- 

From equation (8) we can also obtain an explicit ex- 
pression for the reflection coefficient r in the steady state 
at probe-field resonance. 



r(st) 



^0 



rio72o 



■ 710720 



(9) 



Using equation (9) we find that in the OFF steady state 

the reflection coefficient is ^Qpp/ro « 0.47. Ideally, the 
reflection coefficient could approach the unit, but for the 
parameters corresponding to our qubit this does not hap- 
pen due to the relatively large dephasing rate, 710 = 
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FIG. 6: On-resonant reflected power in the steady state (a) with cross-coupling, and (b) without cross-coupling. 

The two circles indicate the steady states Popp (red) and PoJ] (green) corresponding to our experiment. The decoherence rates 
are 710 = 27i x 14 MHz and 720 = 27r x 27 MHz. 
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FIG. 7: (a) Real part and (b) imaginary part of the reflection coefficient as a function of detuning Ap, at 
different control field amplitudes Qc. 



2rio. Now, when the coupling field is ON and the system 
reached the stationary state Pq^' equation (9) predicts 

that there is almost no reflection (^qn/^o < 0.05) of the 
probe field (green circle in Fig. 6). These results are in 
reasonably good agreement with the numerical simula- 
tions in time domain (see also Fig. 8), even though the 
approximation iytp/Vtcf' ^ 1 is not valid stricto sensu. 
For example, equation (9) predicts that the ON/OFF ra- 
tio i^oFp/^ON is of the order of 100, in agreement with 
the numerical result obtained before. To obtain a higher 
reflection in the OFF state, equation (9) shows that one 
needs to have a lower dephasing rate and to reduce the 
value of the coupling field in the OFF state. If the cou- 
pling field in the OFF state is well below the value of 
the decoherence rates, we have from equation (9) that 



^off/^on ^ (^1 + ^c(ON)/7io720 j . As noticed before, 
increasing l^c(ON) is ultimately limited by anharmonicity. 
Experimentally, lowering the OFF field to zero is limited 
by the leakage of the mixer used to pulse the coupling 
field, while the reduction of dephasing is a general prob- 
lem in the field of superconducting qubits. 

Time- domain transients In Fig. 8 we present the results 
of time-domain numerical simulations of the reflection 
coefficients for a switching event corresponding to the 
one occurring in the experiment presented in Fig. 3. We 
see that, as in the case of the population number, the 
transients tend to stabilize after times of the order of 100 
ns. 

To conclude, our results demonstrate that a super- 
conducting three-level system with large enough anhar- 
monicity and low dephasing can be used for example as 
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FIG. 8: (a) Real part, (b) imaginary part, and (c) 
Smith chart of the normalized reflection coefficient 
in time domain, (d) The corresponding (normalized) 
power reflection. The control field is OFF before the 150 ns 
delay time (red), then switched ON between 150 and 300 ns 
(green), and finally switched OFF again (blue). The dashed 
arrows in (c) indicate the direction of time- evolution. The 
simulations correspond to the experimental data along the 
arrow from Fig. 3. 



the node of a quantum network, redirecting single pho- 
tons along chosen paths (reflecting them in the OFF state 
and allowing them to pass in the ON state). 



Methods 

To probe the qubit we use a weak continuous probe 
tone. The Hamiltonian equation (1) is dynamically ma- 
nipulated by changing the amplitude of the control field, 
resonantly coupling the two higher levels |1) and |2). This 
field is turned on at 150 ns and turned off at 300 ns. For 
measuring the occupation probabilities, we use the fact 
that the three levels reside in a met ast able well . The 
potential barrier that creates this well is then lowered 
by a measurement pulse applied after a measurement 
delay time (see Fig. 2), such that the states above the 
ground state are allowed to tunnel out of the well with 
a certain tunneling probability. The resulting variation 
in the magnetic flux piercing the qubit loop is detected 
by a nearby dc SQUID, which is switched in the run- 
ning state . Therefore from the switching statistics we 
can calculate the total state occupation probability of the 
two excited levels P|i) + ^|2)^^ that have tunneled out. 
The numerical simulations are done by solving the master 
equation with the fourth-order Runge-Kutta method^''. 
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